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The Kennicutt-Schmidt Star Formation Relation at z ^ 2 



o 

Oh, 
< 

O 
(N 

o 
u 

6 



(N 
> 

O 
(N 

O 
rn 

in 

o 
o 



Desika Narayanan^*!, Thomas J. Cox^J, Christopher C. Hayward^, Lars Hernquist^ 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, Ma 02138 
^Observatories of the Carnegie Institution of Washington, 813 Santa Barbara St., Pasadena, Ca, 91101 



MNRAS, 412, 287 



ABSTRACT 



Recent observations of excited CO emission lines from z ^ 2 disc galaxies have 
shed light on the SFR oc p^ relation at high-z via observed Ssfr ~ ^coj=2-i ^^'^ 
SsFR — ^coj=3-2 relations. Here, we describe a novel methodology for utilising these ob- 
servations of high-excitation CO to derive the underlying Schmidt (SFR ex p^) relationship. 
To do this requires an understanding of the potential effects of differential CO excitation with 
SFR. If the most heavily star-forming galaxies have a larger fraction of their gas in highly ex- 
cited CO states than the lower SFR galaxies, then the observed molecular Kennicutt-Schmidt 
index, a, will be less than the underlying SFR oc p^ index, N. Utilising a combination of SPH 
models of galaxy evolution and molecular line radiative transfer, we present the first calcula- 
tions of CO excitation in z ^ 2 disc galaxies with the aim of developing a mapping between 
various observed Esfr — Sqq relationships and the underlying SFR ex p^relation. We find 
that even in relatively luminous z ^ 2 discs, differential excitation does indeed exist, resulting 
in a < iV for highly excited CO lines. This means that an observed (e.g.) Esfr — ^co 1=3-2 
relation does not map linearly to a Ssfr — '^b.2 relation. We utilise our model results to 
provide a mapping from a to N for the range of Schmidt indices N ^ 1 — 2. By com- 
paring to recent observational surveys, we find that the observed Ssfr ~ ^coj=2-i ^^'^ 
SsFR — Scoj=3-2 relations suggest that an underlying SFR oc p^'^ relation describes z ^ 2 
disc galaxies. 

Key words: galaxies: star formation-galaxies:formation-galaxies:high-redshift- 

galaxies:starburst-galaxies:ISM-galaxies:ISM-ISM:molecules-comology:theory 
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1 INTRODUCTION 

Since the original works bv lSchmidJ ( 1 19591) anci lKennicut^ ( Il998bl) 
parameterising star formation rates (SFRs) in galaxies in terms of 
the scaling relation: 



SFR oc pg 



(1) 



there have been considerable efforts by both the Galactic and ex- 
tragalacti c star formation co mmunities to characterise the exponent 
A'^ (e.g. iKeimicutd 1 1998al . and references therein). Constraining 
the Schmidt SFR relation in galaxies is desireable, both for un- 
derstanding the physics of star formation on local scales, as well as 
for giving simulators recipes for modeling processes below typical 
numerical resolution scales. 

Because the volume density is not observable, most observed 
forms of Equation [T] have been in terms of the SFR and gas sur- 



face densitieo As a result of recent pioneering high resolution 
millimetre-wave surveys, two trends have become apparent. 

First, the SFR surface density in galaxies is well correlated 
with the molecular gas surface density, thoug h little relation ex- 
ists between the SFR and HI atomic gas (e.g. Isigiel et aLlllOOa : 
iLerov et alj2008l : ?). Second, the global integrated relationship be- 
tween SFR and CO molecular gas surface densities in observations 
of local galaxies carries an exponent A^ ^ 1.5 (e.g. Sanders et alj 
Il99ll : lKennicuttllT998bl : lOao & Solomon! l2004al llj). A variety of 
theories have been propo sed to explain the Af ^ 1.5 index inferred 
from observations (e.g. ISilldll997l:[Tanll2q0(]|:lElmegreerJ|2002bl : 



iKrumholz & McKee|[2005l : lKrumholz et al.P 2009b^■ most of which 
rely on stars forming on a time scale proportional to the dynamical 
time scale. 
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^ To remain consistent with typical literature nomenclature, we will refer 
to the volumetric form of Equation [T] as the Schmidt relation, and the sur- 
face density form as the Kennicutt-Schmidt (KS) relation. We will reserve 
the index. A', for the volumetric exponent, and a for the surface density 
exponent. 



2 Narayanan, Cox, Hayward & Hernquist 



The recent advent of sensitive (sub)mm-wave interferometers 
has allowed, for the first time, observational constraints on the 
SFR oc p relation in "normal" star-forming (e.g. not starburst- 
ing) galaxies at high redshift via t he detection of CO rotational 
emission li nes iBoucheetal. 2007; lono et alj 20091 : Daddi et al. 
20091120 idllSothwell et al..2009i;iGenzel et alj2010l:lTacconi et'al 



201(J) . However, owing to fact that most mm-wave interferometers 
operate in the 1-3 mm wavelength regime, most detections of CO 
emission lines at high-redshift are of relatively excited lines (e.g. 
CO (J=3-2), as opposed to the ground state transition which is typ- 
ically used as a tracer of H2 molecular gas. 

In principle, if the excitation of CO is relatively invariant in a 
given sample of galaxies, one can simply make an assumption re- 
garding the CO (J=3-2) to (J=l-0) scaling ratio in high-z galaxies, 
and derive an H2 gas mass (or surface density) with the inferred 
CO (J=l-0) intensity. However, if the excitation of CO in galax- 
ies is not constant with increasing SFR, then the exponent in the 
SFR (X CO" relation may not naturally translate to an exponent in 
the SFR (X p^ relation. In other words, a does not always equal N 
when probing high excitation CO lines. 

To see this, consider a sample of galaxies which are forming 
stars at a rate according to SFR oc p^'^. If the CO (J=l-0) line 
faithfully traces the H2 gas mass, then one could expect an ob- 
served relationship SsFR oc S noj^i_o JKrumholz & Thompson! 
I2OO7I ; iNaravanan et al.ll2008cr) . However, if the galaxies with the 
highest SFR have a higher fraction of their CO gas excited into the 
CO J=3 state than the lowest SFR galaxies, one will observe a. flat- 
ter relationship between SFR and Lco3-2 than index a — 1.5. 
The EsFR — Sqoj=3-2 relation will only map linearly to the 
SsFR — 2^2 relation if the excitation of CO is invariant with SFR. 
In practice, this can only happen if the CO gas is thermalised in the 
observed lines (if a given excited state is in LTE). 

This effect of differential molecular excitation on observed 
molecular SFR scaling relations has been observed in the local 
Universe. The SFR-CO (J=l-0) relationship has an index a ~ 
1.5, while the SFR-CO (J=3-2) relationship h a s a flatter index 



0.9 llSandersetal.lll99ll: 



. . lYaoet"aLll2003l : lNai-avananetal.l 

I2OO5I ; llono et alj l2009l : iBavet et alj l2009h . Similarly, while the 
SFR-HCN (J=l-0) relationship is linear in local galaxies, the 
SFR-HCN (J=3-2) index is decidedly sublinear (with index 



0.7; iGao & Solomon 2004a, b; Bu ssmann et al.l |2008| ; 



lOracia-Carpio et alJi2008i ; 1 Juneau et alji2009i') . Observations of in- 
dividual star-forming clumps (which are massive enough to host 
stellar clusters) have been inconclusive regarding whether these 
global trends extend to smaller scales (Wu et al. 2005 , 201Q) . 

In the absence of a more direct tracer of H2 gas than high- 
excitation CO, the potential effects of differential molecular exci- 
tation with SFR need to be quantified in order to derive the un- 
derlying SFR relation in high-redshift galaxies. The few multi- 
line constraints of excitation in high-z galaxies hints that CO may 
be su bthermally excited even in the most luminous z ^ 2 sys- 
tems ('WeiB et al.'2007l; lDannerbauer et alj2009l ; ICarini et al.l2O10l ; 
[Harris et al. 2010). This indicates that applying a uniform mapping 
from (e.g.) CO (3-2) to CO (J=l-0) line intensities will indeed be 
problematic. In this arena, numerical models can offer guidance. 

Our aim in this paper is to calculate the mapping of observed 
SsFR — Sqo relations of excited lines (e.g. CO J=2-l and CO 
J=3-2) to an underlying SFR oc p^ rel ation cont rolling the star 
forma tion. In lNaravanan et alj ( [2009al lbl l2010l) and lHavward et al.l 
1 I2OIOI) , we have developed a merger-driven model for the forma- 
tion of high-redshift ULIRGs which shows reasonable correspon- 
dence with observed SEDs, CO emission properties and number 



counts (C. Hayward et al. in prep.). Here, we utilise the (idealised) 
progenitor disc galaxies of these model mergers to represent the 
star-forming disc s at high- z typically observed in CO emission 
line surveys (e.g. iTacconi et al. 2010: Daddi et al. 2010). We com- 
bine these hydrodynamic simulations of disc galaxies with 3D non- 
LTE molecular line radiative transfer calculations in order to cal- 
culate the full statistical equilibrium excitation properties of the 
molecules. These methods allow us to determine the differential 
excitation of CO of 2 ~ 2 disc galaxies with respect to SFR, and 
derive a mapping of an observed molecular Kennicutt-Schmidt law 
(EsFR — Sgo) to a SFR oc p^ relationship. 



2 NUMERICAL METHODS 

Generally, our goal is to simulate galaxies hydrodynamically that 
serve as reasonable representations of the star forming galaxies at 
z ~ 2 residing on the "main sequenc e" of the SFR-M* relation (e.g. 
not starbursts; iNoeske et al.ll2007 al ibi). We first describe the hydro- 
dynamic methods employed, and then follow with our parameter 
choices which ensure the physical properties of the model galaxies 
are comparable to those observed. 

We simulate the hydrodynamic evolution of the gas phase 
of our model galaxies utilising the fully entropy-conserving A*'- 
body/SPH code GADGET-3 (Springel 2005). The main components 
of the code relevant to this study are the ISM and star formation 
prescriptions. 

The ISM is modeled as multiphase in nature, with cold 
clouds embedded in a hotter, pre ssure-confining ISM (e.g. 
JMcKee & Ostrikedll977l ; ISDringel & H ernquist 2003). This ISM is 
pressurised against runaway fragmentation by supemovae which is 
handled via an effective equati on of state (EOS). For details regard- 
ing this EOS, see Figure 4 of ISpringel et alj ( 120051 ) . Here, we as- 
sume the stiffest equation of state (qeos = 1 in Sprin gel et al.l 
2005). Because the global excitation of CO is to some degree 
dependent on the density structure of the ISM in the galaxy, 
we explore relaxing this assumption with test cases of qeos ~ 
0.75, 0.25. We discuss the magnitude of uncertainty this parame- 
ter causes in the next section. 

Star formation is controlled by a volumetric Schmidt-relation 
such that SFR oc p^ . In order to explore the effects of this index on 
the observed SFR-Lmoi relations, we have run simulations varying 
A'^ between 1, 1.5 and 2. We choose these values as representa- 
tive of the t ypical range of observed Kennicutt-Schmidt indices at 
high-2 (e.g. JBouche et alj|2007l ; JBothwell et ai]|2009l ; JDaddi et al.l 
l20IOl ; lGenzel et al.ll2010l) . The normalisation is anchored such that 
the a volumetric Schmidt relation with index A'^ = 1.5 returns a 
surface- density Kennicutt -Schmidt relation consistent with obser- 
vations JKennicuttll 1 998bl) . The normalisation of the N = 1 and 2 
volumetric relations is forced to match the A'^ = 1 relation at 20 
times the density threshold for star formation. This normalisation 
was chosen such that a disc modeled after Milky Way parameters 
would have a global SFR of ~ 2 MQyr^^ for all three cases. As 
we will discuss later, we explore the effects of varying this normal- 
isation by a factor of 2 in either direction. 

It is worthwhile to consider that we force the stars to 
form according to a volumetric Schmidt relation, though a sur- 
face den s ity K ennic utt-Schmid t relati on is what is observable. 
ISpringell ilOOd) and ICox et alj ( I2OO6I) have shown that a given 
Schmidt (volumetric) relation and Kennicutt-Schmidt (surface den- 
sity) relation have the same exponent in disc simulations very 
similar to these, and we have confirmed this equivalence with 
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the simulations employ e d her e. This similarity is not obvious. 
ISchave & Dalla Vecchial ( 120081) have shown that a linear mapping 
between Schmidt and Kennicutt-Schmidt SFR relations in numer- 
ical simulations is dependent on the choice of equation of state. 
The equation of state we emp loy (particularly the stif f est on e, 
(/EOS = 1) is quite similar to ISchave & Dalla Vecchial J2008r) 's 
"preferred" EOS which reproduces the linear mapping between the 
underlying Schmidt relation and observed Kennicutt-Schmidt rela- 
tion. We therefore refer to the volumetric and surface density rela- 
tions interchangeabljQ 

The simulations are not cosmological. We simulate idealised 
discs in order to maximize spatial resolution (here, the gravitational 
softening length was 100 h^^pc for baryons and 200 h^^pc for 
dark matter). Because the simulations are not cosmological, we ne- 
glect potential gas replenishment from the IGM. However, we are 
less concerned with the temporal evolution of the model galaxy, but 
rather we aim to have the galaxy pass through phases during which 
it has physical parameters comparable to those inferred &X z ^2. 
We thus initialise our simulations with /g=0.8, and allow the disc 
to stabilise for some time. We then consider the snapshots in our 
model discs which have gas fractions between /g=0.4-0.2 as mo- 
tivated both by measurements o f 2: ~ 2 galaxies of this baryonic 
mass (c.f. § [121 lErb et al.ll200^ : lOaddi et alj[2009l ; iTacconi et all 
l2O10f) . as well as the typical steady-state gas fraction of galaxies 
above Mbar > 10^^ in cosmological simulations (e.g. lOave et al.l 
I2OIO !). In a cases where the gas fraction remains above fg > 0.2, 
though the bulk of the gas is below the star formation threshold, 
we arbitrarily cut out snapshots below SFR < 5Mf7^yr ~^ to remain 
consistent with the SFRs of observed galaxies at z ~ 2 JDaddi et al.l 
[20T3). 

W e model discs inside dark matter halos with iHemauisg 
( 11990) density distributions of mass ~ 3 x 10^"^ Mq . The haloes 
are popu l ated w ith discs who are constructed according to the 
IMo et alj ( Il998h formalism. These galaxies are bulgeless and have 
a total baryonic mass of ~ 2 x 10^^ M(Ti, comparable to massiv e 
5zA: galaxies at 2 ~ 2 (e.g. lOaddi et alj2007l ; lTacconi et al.l20ldh . 

To summarise the galaxy evolution modeling, the galaxy snap- 
shots utilised in this study are "selected" according to particular 
criteria in order to represent massive high-z discs. In particular, the 
modeled galaxies have baryonic masses Mbar ~ 2 x IO^^Mq and 
gas fractions /g = 0.2 — 0.4. When analysirig the synthetic SEDs 
from these galaxies, the modeled BzK col or jj are consistent with 
selection as a star-forming BzK galaxy (e.g. lOaddi et alj2005r) . The 
result of this is a galaxy sample which lies on the high-mass end of 
the z ~ 2 SFR-M* relation (Daddi et al. 2007), consistent with the 
typical galaxies ob served for CO emission (e.g. iDaddi et alj20ld : 
iGenzel et alj201(]|) . This selection returns roughly ~30-40 galaxies 
each for the N = 1, 1.5 and 2 model types. 

The simulation snapshots are analysed in post-processing 
with TURTLEBEACH in order to calculate their synthetic CO 
emission properties. TURTLEBEACH is a 3D non-LTE molecular 



^ It is conceivable that in nature these relations ai'c not in fact equivalent. 
However, as this is presently unconstrainable, we are forced to operate un- 
der the assumption that the volumetric and surface-density exponents are 
indeed equivalent in real galaxies. 

^ The sy nthetic colors are calculated with th e dust radiative transfer code 
SUNRISE jjonssonllioo^ Ijonsson et alj|2010l) . We refrain from describing 
the parameters for the dust modeling as reporting the results from these 
calculations is not the primary objective of this paper. The SUNRISE param- 
eters utilised are identical to those in Narayanan et al. (,2010.) , and we refer 
the interested reader to that paper for more details. 



line radiative transfer code which considers both radiative and 
coUisional (de)excitation in c alculating th e excitation c onditions 
dNaravanan et al.l2006ll2008 J) based on the lBernesI ([1979') method. 

The process of calculating the excitation conditions in the 
molecular gas is as follows. First, the level populations (e.g. num- 
ber of molecules at a given excitation level) are guessed at. Based 
on this, each molecular cloud in the galaxy emits CO line photons 
isotropically, with directions determined via Monte Carlo draws. 
These photons are absorbed by molecules in neighbouring clouds. 
Once the mean intensity, J,j, is known across the grid, the coUi- 
sional rateO of CO with H2 (based on the gas densities) are cal- 
culated, and the level populations are updated. At this point, model 
photons are re-emitted based on the new excitation conditions, and 
the whole process is repeated. This process is iterated upon until 
the level populations are converged. 

We assume that half of the cold gas in the GADGET-3 simu- 
lations is in molecular form (as motivated by observations of local 
galaxies; iKeres et al.ll2003h . By assuming that a constant fraction 
of the star-forming gas in the GADGET-3 simulations is molecular, 
we are effectively assuming that the Schmidt law we impose in the 
SPH simulations is a molecular Schmidt law, which we expect is 
reasonable. Within the H2 gas, we ass ume the CO has a uniform 
Galactic abundance of 1.5 x 10"*/H2 dLee et al.lll996h . 

The SPH simulations have resolution of ~100 pc; as such, 
we do not have information regarding the state of the H2 gas be- 
low these scales. Thus, subgrid techniques are necessary to ac- 
count for the density distributions of molecular gas. We assume the 
H2 gas in each cell is bound in giant molecular clouds (GMCs) 
with masses rand omly drawn from the Galactic mass spectrum 
dBlitz et al J 120070 . and density distributions following power-law 
spheres. The clouds follow a Galactic mass-radius relationship (e.g. 
fSolomon et al. 1987 : Rosolowskv 2005. 2007) . Together, these pa- 
rameters set the density distribution of H2 gas in the model galax- 
ies. Because the excitation of molecules is sensitive to the den- 
sity distribution of gas in the galaxies, we explore the effect of 
these parameters on our final results. Observations of GMCs sug- 
est a range of power-law indices, ranging from n = 1 to 2 



Andre et alJl99a : lFuller & Mversll992l : lwalker et aljl990l) . Sim- 
ilarly, the index on the GMC mass spectrum is thought to vary from 
7 sa —1.4 to —2.8 (e.g . lElmegreen 2002a). Tests performed by 



[Narayanan et alj ( l2008a) have shown that within these ranges, the 
CO excitation in galaxies similar to those presented here is not sen- 
sitive to choice of mass-spectrum index or cloud power-law index. 
Nominally, we employ n = 1.5 for the cloud power-law index 
and 7 = —1.8. A more detailed description (including the under- 
lying equations) regard ing our subgrid methods can be found in 
[Narayanan et alj d2008d) . Fi nally, we have benchmarke d our codes 
against literature s tandards ([van Zadelhoff et al J [2002[) . and pub- 
lished the results in [Naravanan et al.[ ( [20060 . 



3 RESULTS AND APPLICATION TO EXISTING 
OBSERVATIONS 

3.1 Main Results 

In Figure[T] we plot the model Esfr — Sgo index, a, as a function 
of CO transition for three model disc galaxies. The disc galaxies 
form stars according to varying Schmidt laws, SFR oc p^, where 



^ The rates come from the Leiden Atomic and Molecular Database 
iSchoier et alj2005l) . 
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Figure 1. Predicted EgpR — SSq index, a, as a function of molecular 
transition for three different galaxy models. The different galaxy mod- 
els vary the relation that controls their SFR such that SFR oc p^ where 
A'^ = 1, 1.5 and 2. The blue shaded region denotes the uncertainties as- 
sociated with limited sample sizes. Generally, the Ssfr — ^noj^i-O 
relation serves as a reasonable proxy for the underlying Schmidt SFR re- 
lation. With higher-lying lines, differential excitation of the molecule with 
SFR becomes an important effect, and the Ssfr — "^qq relation is flat- 
tened from the underlying SFR oc p^ relation. For high critical density 
tracers, a ^ Af. The green star shows the recent CO (J=2-l) data from 
IPaddi et al.l 120101) . and the red star sho ws the CO (J=3-2) measurements 
of 2 ~ 2 galaxies bv lGenzel et al.l 120101) . This data suggests an underlying 
Schmidt index of Af = 1.5. 



A'^ = 1,1.5 and 2. The top curve denotes the N — 2 case, middle 
A*' — 1.5, and bottom TV = 1. The indices, a, are derived by fitting 
log(SFR) versus log(/co) for a random draw of 20 galaxies (corre- 
sponding to typical sample sizes in the current literature) from our 
parent samples of ~ 30 — 40. We do this 1000 times for each model 
and plot the standard deviation in the derived a-indices as the blue 
shaded region. We note that there is an uncertainty of ~ 25% in 
all mean values shown in Figure [T] based on varying initial condi- 
tions. To avoid cluttering the figure, this uncertainty is not denoted 
in Figure[T]itself, though will be discussed later when we provide a 
quantitative mapping between the Schmidt and Kennicutt-Schmidt 
indices. The bulk of the arguments made in this paper are sum- 
marised in this figure. 

First, we see that CO (J=l-0) generally serves a good tracer 
of the H2 molecular gas, and that the Esfr — Ecoj=i-o rela- 
tionship maps reasonably well from the underlying SFR-p^ rela- 
tion. This is b ecause CO (J=l-0) has a relatively low critical density 
( lEvanslll999l) , and most of the molecular gas emits the J=l-0 line. 
While CO (J=l-0) is typically optically thick within GMCs, owing 
to large velocity gradients it is optically thin on galaxy-wide scales. 
Thus so long as GMC properties do not vary strongly from galax 



to galaxy (as loca l measurements suggest: ISolomon et al 



iRosolowskvl I2OO5I: iBlitz & Rosolowskvl I2OO6I: iBlitz et al.l \200lh . 
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1 



and their mass and radius distributions are relatively narrow, then 
an increase in H2 surface density will correspond to an increase in 
the number of GMCs, and a commensurate increase in CO (J=l-0) 
emission. For this reason, the Esfr — Egoj=i-o relation serves 
as a reasonably good proxy for the Schmidt relation. 

The situation is different for higher excitation CO emis- 



sion lines (e.g. CO J=3-2). Higher lying lines have relatively 
high critical densities. For example, the CO (J=3-2) line re- 
quires typical densities of n > lO^cm^^ to excite the line. Be- 
cause of this, these lines are typically su bthe rmal even in the 
most luminous high-redshift galaxies (e.g. 



2003; 



Papadopoulos & Ivison' '2002'; 'Grev e et alj , 

2006; Weifi et al. 2007; Dannerbauer et alj2009l ;ICarilli et a lj20K 
Harris et al.l 120101) , and consequently do not trace the bulk of 
the molecular gas. For higher-lying CO lines, the emission 
line luminosity increases superlinearly with increasing mean 
gas density (owing to the combined effects of increased gas 
mass as well as increased excitation), and the observed SFR- 
Lco" in dex, a, will be less than the underlying Schmidt in- 
dex, A^ jKrumholz & ThomDsonI |2007|; iNaravanan et al.l l2008cl; 
IJuneau et alj2009l) . 

Another way of saying this is that there is a differential ex- 
citation for galaxies with increasing SFR (or mean gas density, 
< n >). The most heavily star- forming (densest) galaxies will 
have a higher (e.g.) CO J=3-2/CO J=l-0 ratio than the lowest SFR 
galaxies. Then, because the SFR-CO (J=l-0) relationship traces the 
underlying Schmidt index, A^, the observed (e.g.) SFR-CO (3-2)" 
relation for higher lying lines will necessarily have a flatter rela- 
tion than the underlying Schmidt relation. This trend will become 
more pronounced as one observes increasingly higher CO (or any 
molecular) transitions with higher critical densities. This is shown 
explicitly in Figure[T] where we see the SFR-Lmoi" index decrease 
as a function of increasing CO transition (or critical density) for all 
model galaxies. Moreover, this was shown to be an observed phe- 
nomen a in the local Universe i n recent surveys bv lBussmann et al.l 
( I2OO8I) and ljuneau sxIaI ( l2009h . 

Figure [T] therefore provides a mapping between the 
Esfr — Ego relation for observed high-lying CO transitions and 
the underlying Schmidt relation. For the purposes of direct appli- 
cation to current high-z data. In Figure|2l we turn Figure[T]into an 
actual mapping between the Esfr — Ego index, a, and the volu- 
metric Schmidt index, A'^. In particular, we focus on the CO (J=3-2) 
transition as it is a relatively commonly observed line at 2 ~ 2. The 
solid line is the mean from Figure[T](e.g. the mean after randomly 
drawing 20 galaxies of our parent sample 1000 times). The blue 
shaded region denotes a 25% range of uncertainty. The uncertainty 
is determined by characterising the dependence of a on the input 
parameters. In particular, within the confines of our "selection cri- 
teria", the input parameters which have the strongest effect on a are 
the equation of state and the Schmidt law normalisation (owing to 
their changing the gas density distribution and the level of thermal- 
isation of the gas). We find the maximum variance in a with these 
parameters is less than 25%. 

We can utilise a combination of recent observations of high- 
z star-forming systems and the information in Figures |2] and [T] 
to infer the Schm i dt relation at high-z. Rec ent investigations by 
iDaddi et al.l ( l201Cl) . lTacconi et alj dloToh and lGenzel et al.l ( |201(I) 
have investigated the Esfr — Ego relation in "normal" star- 
forming systems at 2:= 1-2 whi c h lie on the main sequence of the 
SFR-M* relation. iDaddi et al.l ( l2O10h find a Esf r - Eg n.T^^_i 
index of a ~ 1.31. Similarly, iGenzel et al.l ( 120 101) find a 
Esfr — Egoj=3-2 index for their data set of a « 1.17. Com- 
paring this to Figure [T] we see that these values correspond to an 
underlying Schmidt index of A'^ = 1.5. 

Figures [T] and [2] therefore provide evidence that an Schmidt 
index A'^ ~ 1.5 holds for high-redshift galaxies. It is important 
to note, however, that there is an uncertainty of ~ 25% in these 
models. Observations of numerous CO transitions will help to nar- 



The KS Law at High Redshift 5 




c 

T3 



W 



1.00 



■M 0.10 



0.01 




Observed E « 



■ E i^Q 3_2 " index, a 



^co 1-0 (normalised units) 



Figure 2. The relationship between the SFR oc p index, A'^, and the ob- 



served SsFR - Sf 



index, a. The red dashed line represents the 



mean result, and the blue shaded region denotes a 25% uncertainty (see 
text for details). The red solid region represents the best fit to available 
CO (J=3-2) data from z ~ 2 discs I Genzel et al. 2010), and the thick- 
ness represents the associated uncertainty. The intersection of the ob- 
served Sspp ~ ^co T=3— 2 index, a and the model results suggests that a 
Schmidt index of Af as 1.5 may appropriately describe 2 ~ 2 disc galaxies. 



row the uncertainty in the derived SFR oc p relation by providing 
additional constraints in the N — a space probed in Figure |2] In 
Table[T] we provide the mapping between Schmidt indices, TV and 
molecular Kennicutt-Schmidt indices, a for four CO transitions. 
These numbers provide enough information to create N ~ a plots 
similar to Figure |2] for transitions other than CO (J=3-2), and test 
the concept that a Schmidt index of Af ~ 1.5 describes star forma- 
tion in high- 2: galaxies. 



3.2 Uncertainties and Dependence on Model Physical 
Parameters 

In Figure[3] we show the simulated Espr — Eco,j=i-o plo* f*"" °ur 
model with a Schmidt index of N = 1.5, and denote the galaxies 
which fall within our selection criteria. The trends in Figure [T] can 
depend on how large of a dynamic range the mock observations 
span. 

The trends in Figure [T] depend on the degree of differential 
excitation, which depends on the dynamic range of SFR surface 
densities (or gas surface densities) observed. For example, con- 
sider the case where observations only probe a limited range of 
extremely high SFR galaxies (here, we take this to mean galax- 
ies with gas fraction 0.4 < fg < 0.8). In the high gas frac- 
tion/high SFR regime, the excitation conditions vary extremely 
rapidly. The most gas-rich, densest systems have much of their 
gas thermalised, whereas lower gas fraction galaxies in this range 
(/g Ri 0.4 — 0.6) begin to contain subthermal gas. While it may 
seem counterintuitive that some high gas-fraction galaxies may 
not be fully thermalised in higher CO lines, it is important to re- 
member that the simulations (as do observations) consider global 
emission, from both low and high-density regimes. Subthermal 
CO level pop ulations have been not e d both i n observations o f lo- 
cal galaxies JNaravanan et al.ll2005L boOSal , llono et alj|2009l) . as 
well as in observations of even the densest, most heavily star- 



Figure 3. The model Espr — ^001=1-0 ■'^l^ti'T' fo'' o™ model with un- 
derlying Schmidt index Af = 1.5. The solid line shows a slope of Af = 1.5. 
The ordinate and abscissa are in normalised units. The green circles show 
the galaxies that satisfy our "selection" criteria culls. Including extremely 
gas-rich (red plus signs) or gas-poor systems (where most of the gas is be- 
low the star formation threshold; puiple crosses) may cause the observed 
KS indices to deviate from those in Figure[T] 



forming submillim etre-galaxies at high-redshift JCarilli et al.l201(J : 
iHarris et al.ll2010 t). The rapidly dropping CO excitation (from the 
higher J states to the ground states) at the highest SFRs causes 
the observed molecular KS relation to steepen from the underlying 
Schmidt relation. In the fiducial case of the A^ = 1.5 galaxy, the 
observed EgpR — Eqoj=i-o index can range from ~1.5 to ~2.7 
as we consider increasing numbers of galaxies with gas fractions 
beyond fg > 0.4 in the fits. 

The molecular KS relation at the other end of the spectrum, 
in gas-poor galaxies, is less clear. On one hand, we can expect that 
it would steepen owing to the bulk of the gas in the galaxy be- 
ing below the SF threshold (see, e.g. ' Bigiel et alj|2008 ). However, 
our models don't consider the possib le destruction of molec ular 
gas in low density environments (e.g. iKrumholz et alj|2009ar) . At 
low densities, the neutral ISM is a mixture of HI and H2. This is 
not captured by our models as we are forced by lack of resolution 
to consider the H2 as a fixed fraction of the neutral ISM mass. In 
reality, if the H2 and/or CO mass also drops at low densities, the 
observed KS relation may not become as steep as Figure [3] would 
suggest. 

The trends predicted in Figure [T] are relatively robust within 
the physical parameter range chosen for the bulk of this paper 
(0.2 < /g < 0.4, with a typical dynamic range in SFR of or- 
der ~ 10). As was just shown, large deviations from galaxies of this 
sort may change the observed mapping from a Schmidt relation to a 
KS relation. That said, our assumed range of gas fractions may ac- 
curately represent real 2 '^ 2 galaxies. For example, cosmological 
hydrodynamical simulations suggest that most galaxies with bary- 
onic masses > a few x IO^^Mq (comparable both to those mod- 
eled, as well as those observed) have steady-state gas fractions of 
order ~0.2-0.4 (e.g. Pave et al. 2010). Recent observations appear 
to come to a simila r conclusion . For e xample, CO measurements 
of BzK galaxies by lOaddi et al.l ( 120091) suggest that the expected 
gas fraction for galaxies of the baryonic mass modeled here are 
fg ^ 0.4. Similarly, by inversion of the Kennicutt-Schmidt rela- 



6 Narayanan, Cox, Hayward & Hernquist 



tion. lErb et alj J2006h suggest gas fractions of massive star-forming 
galaxies of order /„ ~ 0.2 — 0.4. Finally, a large CO survey by 
IXacconi et al.l ( 120 lOh suggest that the distribution of gas fractions 
of normal, star-forming galaxies in this mass range between z=\- 
3 is relatively sharply peaked at fg = 0.3 — 0.4 with only rarer 
excursions outside of this range. 

It is additionally worth discussing the potential effects of our 
ISM assumptions on our modeled results. As was discussed in §[2] 
the hydrodynamic galaxy evolution calculations cannot resolve gi- 
ant molecular clouds, the primary origin site of CO emission lines. 
Because of this, we are forced to utilise a subgrid prescription for 
including GMCs in the radiative transfer calculations. In the ab- 
sence of any observational constraints of the nature of the struc- 
ture of the molecular ISM at high-jz, we assume that clouds exist 
as bound spheres, following a Galactic mass-spectrum and Galac- 
tic mass-radius relation with a power-law density gradient (see 
iNaravanan et alj2008bl . for the actual underlying algorithms). 

Our model re sults are dependent on this assumption. As 
was sho wn in .Krumholz & Thompson! ( 120071) and lNaravanan et al.l 
( l2008a) . the degree of differential excitation in galaxies comes 
from the shape of the high-density tail of the molecular gas den- 
sity distribution. While determining the exact dependence of the 
molecular excitation on the density distribution of molecular gas is 
worthy of an independent study, previous model results suggest that 
a range of reasonable assumpti ons result in similar differential ex- 
citation patterns. For example, iKrumhoIz & Thompson! ([20071) as- 
sumed a log-normal dens i ty dist ribution in their model galaxies, 
whereas INaravanan et alj ( l2008c!) utilised a numerical model de- 
rived from the aforementioned cloud mass spectrum, mass-radius 
relation, and power-law density distribution within clouds. Both 
studies found a similar degree of differential excitation in local 
galaxies (see, e.g., the similar predictions for the SFR-HCN 3-2 
relation in local galax ies between the two models in Figure 2 of 
iBussmann et alj|2008 !). In this sense, our model assumptions for 
the structure of the molecular ISM appears reasonable, at least in 
the context of local galaxieo In a similar vein, we note that the 
assumed equation of state in our models can affect the density dis- 
tribution via pressurisation of the ISM. As noted in § 13.11 within 
the ra nge of our modeled EOSs (qeos = 0.25 — 1, see ISpringel! 
I2005I . and §|2]for more details), we find < 25% difference in the 
modeled molecular KS indices. 

It is conceivable, however, that the molecular gas in high- 
redshift galaxies is different in its structure than clouds in the 
Galaxy. For example, observations of local ULIRGs show that in 
high gravitational potentials, tides may cause the molecular ISM to 
exist in a smooth stru cture with a large volume- filling factor, rather 
than in bound clouds dOownes & Solomon!ll998i) . It is not entirely 
clear, in these starburst galaxies, whether the density distribution is 
similar in shape to those in normal discs (e.g. with a high-density 
tail), or whether they have a substantially different gas density dis- 
tribution. Resolved maps of GMCs in local ULIRGs with ALMA 
will help to answer this question. 

Finally, we note that galaxies at high-z may have a larger frac- 
tion of their gas in a dense phase than local galaxi es. This is explic- 
itly acc ounted for in our modeling. As described in lNaravanan et al.l 
( l2008b!) . clouds are randomly drawn from the mass spectrum to fiU 



cells of a given mass until the mass-budget is used. When the den- 
sity in a particular region is high (as informed by the hydrodynamic 
models), more high-mass clouds will statistically be drawn, thus in- 
creasing the dense gas fraction. 

It is important to note, however, that the effects of radiative 
transfer on scales below our cell-sizes are only marginally ac- 
counted for in these models. When a photon is emitted, it sees a 
column density drawn randomly from the distribution of columns 
seen in the cell (Naravanan et al. 2008b), and deposits some inten- 
sity. However, because the velocity distribution is not modeled on 
sub-grid scales, this is a limiting assumption. Some of this radiation 
may actually emerge from the cell owing to large velocity gradients 
within a cell which may keep photons from being trapped. Efforts 
are underway to explicitly account for this effect, though it is a task 
well outside the scope of this study. 



4 DISCUSSION 

The principle result of this study is that, due to subther- 
mal excitation in high-lying CO lines in high-z galaxies, ob- 



served Ss 



■^co 



relations do not necessarily map linearly 



^ It is additionally worth noting that [Narayanan et alj <2008cl) varied the 
range of indices of the Galactic mass-spectrum and power-law density gra- 
dients within the range of observational constraints, and found little differ- 
ence in the final differential excitation in their modeled local galaxies. 



to EsFR — ^§2 relations. T his is part of a broader theore tical 
framework first develope d bv lKrumholz & Thompson ! (lOfft) and 
INaravanan et alj ( l2008q) which posits that the underlying den- 
sity distribution of galaxies is crucial in mapping the underly- 
ing Schmidt relation to observed molecular Kennicutt-Schmidt 
relations. In this paper, we have expanded upon these studies 
by exploring the relationship between Schmidt and Kennicutt- 
Schmidt indices in models which aim to serve as reasonable rep- 
resentations of the 2 ~ 2 star-forming disc gal axies being uncov- 
ered in sensitive optical/NIR observations (e.g. lOaddi et al.ll2005! : 
iForster Schreiber et al.l2009l) . Indeed, the models studied here have 
been shown in previous publications to satisfy the star-forming 
BzK colour-selection criteria (Narayanan et al. 2009b), as well as 
serve as reasonable progenitors for luminous , merger-driven z ~2 
submillimetre and 24 yim- selected galaxies JHavward et alj|2010l : 
INaravanan et alJuOlOf) . We have additionally explored the depen- 
dence of the observed molecular Kennicutt-Schmidt index on vary- 
ing underlying Schmidt relations. Our models, combined with the 
gen eral theoretical framewor k of iKrumholz & Thompson! <2007r) 
and INaravanan et al.l ( l2008c !) suggest that observed z ^2 discs 
are subject to differential CO excitation with respect to SFR, and 
that the observed Esfr - Egoj=2-iand Esfr - Egoj=3„2 re- 
lations may map to an underlying Schmidt index of A'^ = 1.5 con- 
trolling the SFR. 

An important verifying aspect to these models is that they are 
able to explain the multitude of observed SFR-Lmoi relations in 
the local Universe. For example, a linear relationship has been ob- 
served between SFR and HCN (J=I-0) in local galaxies, a trend 
which has often been interpreted as evidence that a linear dense gas 
Schmidt relation controlled the SFR. Our models suggest the linear 
SFR-HCN (J=I-0) relation is in reality a combined effect of an un- 
derlying S FR PC p^ index of N = 1.5 and differential excitatio n 
in HCN tKrumholz & Thompson! 1 2007I ; INaravanan et alj|2008c!). 
This v iew has been observationally confirmed by iBussmann et al.l 
( 12008!) . who showed that local galaxies exhibit an sublinear SFR- 
HCN (J=3-2) relation, and thus follow the trend of decreasing SFR- 
HCN° index with increasing transition number characteristic of 
these models (e.g. Figure [TJ see also Figure 7 of INaravanan et al.l 
l2008q) . Similarly, this model satisfies the multi-line constraints of 
local CO observations. The SFR-CO (J=l-0) relation in local galax- 
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Table 1. This Table provides the mapping between observed Kennicutt-Schmidt molecular surface-density indices, a, and underlying volumetric Schmidt 
indices, N . We provide the mapping for CO transitions J=l-0 through J=4-3. The numbers contained here constitute the information necessary to re-create a 
plot like Figure|2]for any of the four modeled CO transitions and will aid interpretation of future observations of varying molecular transitions. The "errors" 
denote a 25% uncertainty level which encompasses variations in the final solution upon changing initial conditions in our simulation (primarily the equation 
of state and normalisation of the SFR relation). 

Schmidt Index Af SFR-CO (J=l-0)" SFR-CO (J=2-l)« SFR-CO (J=3-2)" SFR-CO (1=4-3)'^ 



1 

1.5 

2 



0.96±0.24 
1.47 ±0.37 
2.13 ±0.53 



0.89±0.22 
1.36±0.34 
1.95±0.48 



0.76±0.19 
1.08±0.27 
1.45±0.36 



0.49±0.12 
0.77±0.19 
0.95±0.24 



ies appears to have an index ranging from ~ 1.3 — 1.5. At higher- 
lying transitions (e.g. CO J=3-2), the i ndex drops to ~ 0.9, in ac- 
cordance with theoretical predictions dlono et al.ll2009h . We note 
that at lower bolometric luminosities, even the SFR-CO (J=l-0) 
relation may be subject to differential excitation and serve as a rel- 
atively poor tracer of the underlying Schmidt relation. 

Finally, with an eye toward ALMA, we comment on the role of 
spatial resolution in observational determinations of the KS relation 
at high-z. The exact mapping between the observed Esfr — Sgo 
index, a, and the SFR oc p index, A'^, depends on the level of ther- 
malisation of the gas within the beam. It is dependent (to first or- 
der) on the mean gas density. Higher resolution observations which 
probe just the nucleus of the galaxy will probe higher mean den- 
sities, and allow even higher-lying CO emission lines to directly 
trace the underlying Schmidt SFR relation. Indeed, some very ten- 
tative observat ional evidence for this t rend in local galaxies has 
been shown bv lNaravanan et alj ( l200 8a). The trends shown in Fig- 
ures [T] and |2] are for the central 6 kpc of the galaxy, comparable 
to the typical resolu tion of current interferometric observations of 
2 ~ 2 galaxies (e.g. iTacconi et alj|2010i" ). Our models suggest that 
observations of the central ~2 kpc will probe sufficiently dense gas 
that the observed Esfr — Egoj=3-2 index, a, will trace the un- 
derlying SFR oc p'^ index, N. 



5 SUMMARY 

Current facilities demand that observations of the molecular 
Kennicutt-Schmidt relation at 2 ~ 2 probe highly excited CO lines 
(e.g. CO J=3-2). However it is not clear how exactly to map ob- 
served Esfr — Egoj=3-2 relations to underlying SFR tx p^ re- 
lations: differential excitation of CO with SFR may make interpre- 
tation difficult. 

In order to aid in the interpretation of observed molecular 
Kennicutt-Schmidt relations, we have calculated the first models 
of CO excitation for star-forming disc galaxies at z ~ 2. Our main 
results are: 

(i) Due to differential excitation of CO with SFR in 2: ~ 2 disc 
galaxies, global observations of a (e.g.) Esfr — Egoj=3-2 rela- 
tion will result in a flatter index, q, than the underlying Schmidt 
(SFR oc p^) index, TV. This trend exists for all lines above CO 
(J=l-0), though the disparity between a and A^ grows with increas- 
ingly high CO transitions owing to increasing critical densities in 
the line. 

(ii) We present a mapping from observed Esfr — Egg indices, 
Q, to underlying Schmidt indices, A'^ for global observations of 
2; ~ 2 discs. Combining these model results with the observed 
(nearly) linear relationship between Esfr and Ecoj=3-2 suggests 
that a relation SFR oc p^'^ may control the SFR in z ~ 2 galaxies. 



ACKNOWLEDGEMENTS 

We thank Shane Bussmann, Neal Evans, Reinhard Genzel, Patrik 
Jonsson, Dusan Keres, Mark Krumholz, Charlie Lada, Kai Noeske, 
Alice Shapley, Amiel Sternberg and Linda Tacconi for enjoyable 
conversations as the ideas for this study were developed. The simu- 
lations in this paper were run on the Odyssey cluster, supported by 
the Harvard FAS Research Computing Group. 



References 

Andre, R, Ward-Thompson, D., & Motte, F 1996, A&A, 314, 625 
Andreani, R, Cimatti, A., Loinard, L., & Rottgering, H. 2000, 

A&A, 354, LI 
Bayet, E., Gerin, M., Phillips, T. G., & Contursi, A. 2009, MN- 

RAS, 399, 264 
Bemes, C. 1979, A&A, 73, 67 
Bigiel, F, Leroy, A., Walter, F, Brinks, E., de Blok, W. J. G., 

Madore, B., & Thomley, M. D. 2008, AJ, 136, 2846 
Blitz, L., Fukui, Y., Kawamura, A., Leroy, A., Mizuno, N., 

& Rosolowsky, E. 2007, in Protostars and Planets V, ed. 

B. Reipurth, D. Jewitt, & K. Keil, 81-96 
Blitz, L. & Rosolowsky, E. 2006, ApJ, 650, 933 

Bothwell, M. S. et al. 2009, MNRAS in press: ArXiv/0912.1598 
Bouche, N. et al. 2007, ApJ, 671, 303 

Bussmann, R. S., Narayanan, D., Shirley, Y L., Juneau, S., Wu, 
J., Solomon, P. M., Vanden Bout, P. A., Moustakas, J., & Walker, 

C. K. 2008, ApJ, 681, L73 

Carilli, C. L. et al. 2010, ApJ Submitted: ArXiv/1002.3838 
Cox, T. J., Jonsson, P., Primack, J. R., & Somerville, R. S. 2006, 

MNRAS, 373, 1013 
Daddi, E., Elbaz, D., Walter, F, Boumaud, F, Salmi, F, Carilli, 

C, Dannerbauer, H., Dickinson, M., Monaco, P., & Riechers, D. 

2010, ArXiv/1003.3889 
Daddi, E. et al. 2005, ApJ, 631, L13 
— . 2007, ApJ, 670, 156 
— . 2009, ArXiv/09 11.2776 
Dannerbauer, H., Daddi, E., Riechers, D. A., Walter, F, Carilli, 

C. L., Dickinson, M., Elbaz, D., & Morrison, G. E. 2009, ApJ, 

698, L178 
Dave, R., Finlator, K., Oppenheimer, B. D., Fardal, M., Katz, N., 

Keres, D., & Weinberg, D. H. 2010, MNRAS, 404, 1355 
Downes, D. & Solomon, P M. 1998, ApJ, 507, 615 
Elmegreen, B. G. 2002a, ApJ, 564, 773 
— . 2002b, ApJ, 577, 206 
Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., Reddy, 

N. A., & Adelberger, K. L. 2006, ApJ, 647, 128 
Evans, II, N. J. 1999, ARA&A, 37, 31 1 
Forster Schreiber, N. M. et al. 2009, ApJ, 706, 1364 



8 Narayanan, Cox, Hayward & Hemquist 



Fuller, G. A. & Myers, P. C. 1992, ApJ, 384, 523 

Gao, Y. & Solomon, P. M. 2004a, ApJS, 152, 63 

— . 2004b, ApJ, 606, 271 

Genzel, R. et al. 2010, MNRAS accepted; arXiv/1003.5180 

Gracia-Carpio, J., Garcia-Burillo, S., Planesas, P., Fuente, A., & 

Usero, A. 2008, A&A, 479, 703 
Greve, T. R., Ivison, R. J., & Papadopoulos, P P 2003, ApJ, 599, 

839 
Hainline, L. J., Blain, A. W., Greve, T. R., Chapman, S. C, Small, 

I., & Ivison, R. J. 2006, ApJ, 650, 614 
Harris, A. I., Baker, A. J., Zonak, S. G., Sharon, C. E., Genzel, R., 

Ranch, K., Watts, G., & Creager, R. 2010, ArXiv e-prints 
Hayward, C. C., Narayanan, D., Jonsson, P., Cox, T. J., Keres, D., 

Hopkins, P. F, & Hernquist, L. 2010, Conference Proceedings 

for UP2010: Have Observations Revealed a Variable Upper End 

of the Initial Mass Function? Treyer, Lee, Seibert, Wyder, Neil 

eds. arXiv/1008.4584 
Hernquist, L. 1990, ApJ, 356, 359 
lono, D., Wilson, C. D., Yun, M. S., Baker, A. J., Petitpas, G. R., 

Peck, A. B., Krips, M., Cox, T. J., Matsushita, S., Mihos, J. C, 

& Pihlstrom, Y 2009, ApJ, 695, 1537 
Jonsson, P 2006, MNRAS, 372, 2 

Jonsson, P, Groves, B. A., & Cox, T. J. 2010, MNRAS, 186 
Juneau, S., Narayanan, D. T, Moustakas, J., Shirley, Y. L., Buss- 

mann, R. S., Kennicutt, R. C, & Vanden Bout, P A. 2009, ApJ, 

707, 1217 
Kennicutt, Jn, R. C. 1998a, ARA&A, 36, 189 
— . 1998b, ApJ, 498, 541 

Keres, D., Yun, M. S., & Young, J. S. 2003, ApJ, 582, 659 
Krumholz, M. R. & McKee, C. R 2005, ApJ, 630, 250 
Krumholz, M. R., McKee, C. R, & Tumlinson, J. 2009a, ApJ, 693, 

216 
— . 2009b, ApJ, 699, 850 

Krumholz, M. R. & Thompson, T. A. 2007, ApJ, 669, 289 
Lee, H., Bettens, R. P A., & Herbst, E. 1996, A&AS, 119, 111 
Leroy, A. K., Walter, F, Brinks, E., Bigiel, F, de Blok, W. J. G., 

Madore, B., & Thornley, M. D. 2008, AJ, 136, 2782 
McKee, C. F & Ostriker, J. R 1977, ApJ, 218, 148 
Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Narayanan, D., Cox, T. J., Hayward, C. C, Younger, J. D., & 

Hernquist, L. 2009a, MNRAS, 400, 1919 
Narayanan, D., Cox, T. J., & Hernquist, L. 2008a, ApJ, 681, L77 
Narayanan, D., Cox, T. J., Kelly, B., Dave, R., Hernquist, L., Di 

Matteo, T, Hopkins, P. F, Kulesa, C, Robertson, B., & Walker, 

C. K. 2008b, ApJS, 176,331 
Narayanan, D., Cox, T. J., Shirley, Y, Dave, R., Hernquist, L., & 

Walker, C. K. 2008c, ApJ, 684, 996 
Narayanan, D., Dey, A., Hayward, C, Cox, T. J., Bussmann, 

R. S., Brodwin, M., Jonsson, P., Hopkins, P., Groves, B., 

Younger, J. D., & Hernquist, L. 2009b, MNRAS accepted: 

arXiv/0910.2234 
Narayanan, D., Groppi, C. E., Kulesa, C. A., & Walker, C. K. 

2005, ApJ, 630, 269 
Narayanan, D., Hayward, C. C, Cox, T. J., Hernquist, L., Jonsson, 

R, Younger, J. D., & Groves, B. 2010, MNRAS, 401, 1613 
Narayanan, D., Kulesa, C. A., Boss, A., & Walker, C. K. 2006, 

ApJ, 647, 1426 
Narayanan, D. et al. 2008d, ApJS, 176, 331 
Noeske, K. G. et al. 2007a, ApJ, 660, L47 
— . 2007b, ApJ, 660, L43 

Papadopoulos, R R & Ivison, R. J. 2002, ApJ, 564, L9 
Rosolowsky, E. 2005, PASR 117, 1403 



— . 2007, ApJ, 654, 240 

Sanders, D. B., Scoville, N. Z., & Soifer, B. T. 1991, ApJ, 370, 
158 

Schaye, J. & Dalla Vecchia, C. 2008, MNRAS, 383, 1210 

Schmidt, M. 1959, ApJ, 129, 243 

Schoier, F L., van der Tak, F F S., van Dishoeck, E. F, & Black, 
J. H. 2005, A&A, 432, 369 

Silk, J. 1997, ApJ, 481, 703 

Solomon, R M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 
319,730 

Springel, V. 2000, MNRAS, 312, 859 

— . 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T, & Hernquist, L. 2005, MNRAS, 361, 
776 

Springel, V. & Hernquist, L. 2003, MNRAS, 339, 289 

Tacconi, L. J. et al. 2010, Nature, 463, 781 

Tan, J. C. 2000, ApJ, 536, 173 

van Zadelhoff, G., Dullemond, C. P., van der Tak, F F S., Yates, 
J. A., Doty, S. D., Ossenkopf, V., Hogerheijde, M. R., Juvela, 
M., Wiesemeyer, H., & Schoier, F L. 2002, A&A, 395, 373 

Walker, C. K., Adams, F C, & Lada, C. J. 1990, ApJ, 349, 515 

WeiB, A., Downes, D., Walter, F, & Henkel, C. 2007, in Astro- 
nomical Society of the Pacific Conference Series, Vol. 375, From 
Z-Machines to ALMA: (Sub)Millimeter Spectroscopy of Galax- 
ies, ed. A. J. Baker, J. Glenn, A. I. Harris, J. G. Mangum, & M. S. 
Yun, 25-1- 

Wu, J., Evans, N., Shirley, Y, & Knez, C. 2010, ArXiv e-prints 

Wu, J., Evans, II, N. J., Gao, Y, Solomon, R M., Shirley, Y L., & 
Vanden Bout, R A. 2005, ApJ, 635, L173 

Yao, L., Seaquist, E. R., Kuno, N., & Dunne, L. 2003, ApJ, 588, 
771 



1000 F 







Pi 

00 




lO'i 
Stellar Mass (M®) 



